The Two-Dimensional Bose-Einstein Condensate 
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We study the Hartree-Fock-Bogoliubov mean-field theory as applied to a two-dimensional finite 
trapped Bose gas at low temperatures and find that, in the Hartree-Fock approximation, the system 
can be described either with or without the presence of a condensate; this is true in the thermody- 
namic limit as well. Of the two solutions, the one that includes a condensate has a lower free energy 
at all temperatures. However, the Hartree-Fock scheme neglects the presence of phonons within the 
system, and when we allow for the possibility of phonons we are unable to find condensed solutions; 
the uncondensed solutions, on the other hand, are valid also in the latter, more general scheme. Our 
results confirm that low-energy phonons destabilize the two-dimensional condensate. 



I. INTRODUCTION 

In a recent paper 0] , one of the groups that pioneered 
the formation and detection of Bose-Einstein condensa- 
tion (BEC) in harmonically trapped atomic gases |2[ 0, ||] 
reports the creation of (pseudo-) two-dimensional con- 
densates. These have been produced by taking a three- 
dimensional condensate of 23 Na atoms and carrying out 
two independent processes on it: 1) Initially, one of the 
confining frequencies (that in the z direction, uj z , in order 
to minimize the effects of gravity) is increased until the 
condensate radius in that dimension is smaller than the 
healing length associated with the interaction between 
atoms (taken here to be repulsive and parametrized by 
the two-body scattering length a). This is not sufficient 
to reduce the dimensionality, however, since the atoms, 
each of which has mass to, will literally squeeze into the 
third dimension if there are more than 
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of them in the trap. 2) Consequently, the number of 
atoms in the condensate must be reduced; this is achieved 
by exposing the condensate to a thermal beam. The 
reduction in effective dimensionality becomes apparent 
when the aspect ratio of the condensate, which is in- 
dependent of N in 3D, starts to change as the number 
of atoms is gradually reduced. The condensates thus 
produced have a number of atoms that ranges between 
10 4 and 10 5 . 

This constitutes an important experimental contribu- 
tion to the long-standing debate about the existence of 
BEC in two dimensions. It is a well-known fact (and 
a standard textbook exercise) that BEC cannot happen 
in a 2D homogeneous ideal gas; a rigorous mathematical 
theorem || extends this result to the case where there are 
interactions between the bosons. When the system is in a 
harmonic trap, on the other hand, BEC can occur in two 
dimensions || below a temperature k B T c — hto yj '6N / n 2 , 
but the theorem is once again valid when interactions 
between the bosons arc considered 0: while there is a 
BEC in the 2D system, it occurs at T = 0. 



The preceding discussion is valid only in the thermo- 
dynamic limit, which in the particular case of an isotrop- 
ically trapped system consists of making N — > oo and 
lu — > in such a way that Nto 2 remains finite |Q |). 
The question remains whether a phenomenon resembling 
BEC — that is, the accumulation of a macroscopic num- 
ber of particles in a single quantum state — occurs or not 
when the system consists of a finite number of particles 
confined by a trap of finite frequency, as is certainly the 
case in experimental situations. If there is such a phe- 
nomenon, one would like to know more about the process 
by which this "condensate" is destabilized at finite tem- 
peratures as N grows. 

Some authors @, 10 have considered, in the finite ho- 
mogeneous 2D case, the possibility of a BEC. A similar 
analysis (ll) was considered for a quasi-2D trapped gas; 
the latter reference finds that the phase fluctuations in 
the condensate vary with temperature and particle num- 
ber as 



((<5</>) 2 ) ocTlogTV 



(2) 



which diverges for finite temperature as N — > oo, as one 
would expect. For finite N, on the other hand, the fluc- 
tuations are tempered at very low temperature; since 
the coherence length, though finite, is still larger than 
the characteristic length imposed on the system by the 
trap ||l2f or by walls, one can speak of a "quasiconden- 
sate." It is this quasicondensation that we wish to study 
in this paper. Reference Jig] reports the observation of 
a quasicondensate in a homogeneous system of atomic 
hydrogen adsorbed on 4 He; at this point we cannot tell 
if the condensates reported in Ref. Q are actually qua- 
sicondensates. 

One can approach the study of the 2D Bose gas by em- 
ploying mean-field theory, which in this context refers to 
the Hartree-Fock-Bogoliubov (HFB) equations and var- 
ious simplifications thereof that we will review in quan- 
titative detail and classify in the present work. The 
HFB theory has been remarkably successful in treating 
the 3D case; though we give a few representative refer- 
ences below, we refer the reader to Ref. 14 , which in- 



cludes a comprehensive review and an exhaustive list of 
references, and will henceforth concentrate on the work 
that has been carried out in 2D. 
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The HFB theory assumes from the outset that the sys- 
tem is partially condensed, and proposes separate equa- 
tions to describe the condensate and the uncondensed 
(thermal) component. The condensate is described by 
a macroscopic wavefunction that obeys a generalized 
Gross-Pitaevskh (GP) equation Jl5|, [Tq] , while the non- 
condensate consists of a superposition of Bogoliubov 
quasiparticle and "quasihole" excitations weighted by a 
Bose distribution. The expression for the noncondensate 
can be simplified by neglecting the quasihole excitations 



, by performing a semi- 
fl , or by combining both 



(the Hartree-Fock scheme) jn| 
classical WKB approximation 

of these [^0], ^l], |22| . Each one of the above schemes can 
be further simplified, when the thermodynamic limit is 
approached, by neglecting the kinetic energy of the con- 
densate: this is the Thomas- Fermi limit [^3[ ^4|. Finally 
one can neglect the interactions between thermal atoms, 
arriving at the semi-ideal model [^5[ ^6) . 

The semi-ideal model has been implemented in 2D J27], 
p8| , but has been found to yield unphysical results as 
the interaction strength becomes sizable. The full-blown 
HFB model does not appear to be much more success- 
ful, at least when considered in the semiclassical limit. 
In a previous paper [ p9[ one of us found that, below a 
certain temperature, the introduction of a condensate in 
the Thomas-Fermi limit (corresponding to the thermo- 
dynamic limit) renders the HFB equations incompatible, 
with the noncondensate density becoming infinite at ev- 
ery point in space. The singularity occurs at the low 
end of the energy spectrum, indicating that the conden- 
sate is being destabilized by long- wavelength phonons; 
this interpretation in terms of phase fluctuations had al- 
ready been proposed for the homogeneous system . 
In this work we also report our inability to find self- 
consistent semiclassical solutions to the HFB model when 
a finite trapped system is considered. 

Moreover, it was recently discovered |3l], [52| that 
it is possible, in the 2D case exclusively, to solve the 
HFB equations semiclassically at any temperature with- 
out even having to invoke the presence of a condensate 
(thus obtaining what we will call an "uncondensed" so- 
lution). In other words, it is possible to simply cross out 
the condensate component and solve for the system to a 
temperature close to T = Q. The solution thus obtained 
shows an accumulation of atoms at the center of the trap 
and yields a bulge in density similar to that caused by 
the presence of a condensate, even though no state is 
macroscopically occupied. 

These results appear to reinforce the conclusion that 
BEC cannot happen in the two-dimensional trapped sys- 
tem. However, we are still confronted with the experi- 
mental results described above. Furthermore, two differ- 
ent Monte Carlo simulations (3^, |5| show significant con- 
centrations of particles in the lowest energy state for fi- 
nite N, though it must be said that this method provides 
little information about the types of excitations that con- 
tribute to the disappearance of the condensate, and that 
it is difficult to carry out such simulations on very large 



systems. 

The HFB model cannot be discounted at this stage 
either: when we restrict ourselves to the Hartree-Fock 
approximation, it is possible to find self-consistent solu- 
tions (henceforth referred to as "condensed") involving 
a condensate, both for finite systems and in the ther- 
modynamic limit, when using either the discrete set of 
equations j35| or the WKB approximation |36 to treat 



the noncondensate. The Hartree-Fock approximation ne- 
glects phononlike excitations, so it is not surprising that 
it yields solutions. However, one may be able to justify 
its usefulness. 

It is known that, in the infinite homogeneous sys- 
tem, infrared singularities occur but are renormalized 
by interactions, providing, in essence, a cutoff at a low 
wavenumber ko 



VnmU |37|, |38|, |3j, where n is the 
density, m the particle mass, and U the effective inter- 
action strength. Indeed, it is possible to estimate the 
Berezinskii-Kosterlitz-Thouless (BKT) transition tem- 
perature by simply cutting off the ideal-gas density ex- 
pression at this fco p8| . Presumably a similar situation 
occurs in the trapped case. The Hartree-Fock approach 
provides a convergent theory by cutting off the singular- 
ities at a wavenumber similar to that of more rigorous 
theories. Whether such an approach gives a reasonable 
estimate of the BKT transition temperature, the super- 
fluid density, or a quasicondensate density for the finite 
system will need to await a more rigorous theoretical ap- 
proach to the interacting 2D trapped gas f4fj[| . 

Given the above limitations, we analyze the char- 
acter of the BEC in the 2D trapped system by solv- 
ing the coupled equations of the theory. We find 
that, in the Hartree-Fock scheme, it is possible to find 
both condensed and uncondensed solutions for the two- 
dimensional equations. We also calculate the free en- 
ergy corresponding to each one and find that the con- 
densed solution has a lower free energy at all tempera- 
tures, which appears to imply that, at least at this level of 
approximation, the uncondensed solution is unphysical or 
metastable. The condensed solution will be "preferred" 
over the uncondensed one, and we assume that the solu- 
tion represents an approximation to a quasi-condensate. 

It is evident from our discussion that our approach to 
BEC in 2D trapped systems is a preliminary one and 
that a further analysis that takes fully into account the 
BKT transition is necessary. A start to answering this 
need has been presented in Ref. Q, and we intend to 
return to this problem ourselves in the future. 



II. THE MODEL 

A. The Hartree-Fock-Bogoliubov equations 

Throughout this paper we use a dimensionless system 
of units in which all len gths are scaled by the oscilla- 
tor length xo = ^h/muj and all energies are expressed 
in terms of the one-dimensional ground-state energy of 
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the oscillator, hu/2. Dimensionless variables will in gen- 
eral carry a tilde: for example, the total density n be- 
comes h = XqU and the chemical potential ft = fi/hhuj. 

The HFB equations Jl|, |l|, || result from assum- 
ing that 1) the (repulsive) interactions between atoms 
consist exclusively of two-body low-energy collisions that 
can be described by a delta- function pseudopotential 
of strength g (related to its dimensionless counterpart 7 
through g = ^Tluixqj), that 2) the many-body field op- 
erator ^ can be decomposed via 



* = (vp) + ip = $ + -0, 



(3) 



where the ensemble average $ is a real macroscopic wave- 
function that describes the condensate (reflecting the 
imposition of macroscopic long-range order on the con- 
densed system), and that 3) products of noncondensate 
operators can be simplified using a finite-temperature 
version of Wick's theorem E3]. If we insert (||) into the 
many-body grand-canonical Hamiltonian 



H = d 2 x^(A 



(4) 



where A = — V 2 +x 2 — ft, and neglect anomalous averages 



via the Popov approximation J2 1[| . we obtain an expres- 
sion that can be diagonalized and yields an infinite set of 
coupled differential equations. 

On one hand, the macroscopic wavefunction mentioned 
above is the square root of the dimensionless condensate 
density n and obeys the generalized Gross-Pitaevskii 
equation, 



A$ + 7(710 + 2n')$ = 0, 



(5) 



where n' = n — ho is the noncondensate density. The 
factor of 2 in (|J) and hereafter is a consequence of the 
direct (Hartree) and exchange (Fock) terms being iden- 
tical, which follows from the fact that the delta-function 
interaction that we are considering has zero range J42| . 
(The term involving the condensate does not include 
that factor; this is a consequence of the restricted grand- 
canonical ensemble that we are using. See the end of the 
next section.) 

The noncondensate, in turn, is described by an infinite 



number of pairs of functions that obey 41 



A + 2771 —jfiQ 
—jfio A + 2jn 

and which generate the noncondensate density via 

«'( x ) = ((N 2 + \ v o\)fj + N 2 )> 



(6) 



(7) 



where we have introduced the Bose-Einstein distribution 
factor fj = (e*^/* — l) 1 , which appears when the free 
energy of the system is minimized, and the dimensionless 
temperature t = k B T/^huj. The last term of (0) describes 
the zero-temperature depletion of the condensate, which 



accounts for less than 1% of the particles and is therefore 
negligible [|T). This self-consistent set of equations is 
closed, and the chemical potential found, by imposing 
that the total number of particles remain fixed: 



N = N, 



+ N' = J d 2 xh(x) = J d 2 x ( 



h + n') . (8) 



At temperatures such that k B T 3> Tiuj, one can use 
the semiclassical (WKB) approximation |19j that results 



from taking Uj 



i(x)e^ and Vj w u(x)e^. The 



phase common to both defines a quasiparticle momen- 
tum through k — V0; u, v, and n vary sufficiently slowly 
with x that their spatial derivatives can be neglected . 
The distribution factor is now fj m /(x, n), infinite sums 
are transformed into momentum integrals — which in two 
dimensions can be solved in closed form |2^| — and the 
equations become algebraic, yielding the Bogoliubov en- 
ergy spectrum |21| 

Z hfb (k, x) = aJ (k 2 + x 2 + 2jh - ft) 2 - 7 2 7t j q (9) 

and the following integral expression for the nonconden- 
sate density: 



t 

47T 

t 
47T 



((x 2 +2=fii-ji) 2 -y 2 hl) 1 / 2 /t cs _ 1 

log(l - e - { \ / ( x2 + 2 ^-^ 2 -"> 2A ?>)/ t ) .(IQ) 



At high enough temperatures there is no condensate, and 
the density on the left-hand side of ( |1C| ) is just the total 
density. We thus have to solve a single self-consistent 
equation, 

n'(x) nfx) = -4- log (1 - e ^ 2 + 2 ^-^/ t ) ; (11) 

47T 

the chemical potential is once again calculated by requir- 
ing N to be fixed. Equation (11) has been found by the 
authors of Ref. |n| to be soluble at all temperatures, 
a result that we confirm in the present work. Refer- 
ence had mistakenly concluded that it was impossible 
to solve (11) below a certain temperature. 

The Hartree-Fock approximation |l7], ^l], |2^| amounts 
to neglecting the Vj in (|6|) and, when combined with the 
semiclassical treatment, results in the disappearance of 
the last term within the square root of Eq. (||). The 
energy spectrum now becomes 

e HF = e = + x 2 + 2jh — ft (12) 

and the noncondensate density turns into 

ft'fx) = -— log (1 - e -(- 2 +27S-A)/t) j ( 13 ) 
47T V ' 

an expression that differs from ([ll]) in that the left- 
hand side corresponds only to the noncondensate density. 
The three-dimensional version of this equation, coupled 
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and show that the latter reduces to 



2 4 6 

x = r/x 
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FIG. 1: Condensate and noncondensate density profiles of 
a two-dimensional gas with N = 10 4 and 7 = 0.1 at T = 
0.7T C , where T c is the ideal-gas transition temperature. The 
Goldman-Silvera-Leggett model (dashed lines) treats the con- 
densate in the Thomas- Fermi limit (|l4|), while the Hartree- 
Fock model (full lines) uses the full Gross-Pitaevskii equa- 
tion (^|). Both models use Eq. (|l^) to describe the noncon- 
densate. 



with (|J) and (|J), has been frequently used in the litera- 
ture to study the three-dimensional gas. Reference [p2| , 
for example, exhibits a detailed comparison of its pre- 
dictions to those of Monte Carlo simulations and hnds 
excellent agreement between the two. 

In two dimensions, the Hartree-Fock equations have 
been solved without resorting to the WKB approxima- 
tion [^H . The authors of this reference succeeded in find- 
ing self-consistent density profiles and used them to study 
the temperature dependence of the condensate fraction. 
Our results |36|, which do take advantage of the semi- 
classical approximation, agree quite well with theirs (see 
Fig.@). 



B. The Thermodynamic Limit 

When N is large we can neglect the kinetic energy 
of the system, which in 2D can be shown to be of or- 
der XlN p9fl , and obtain the Thomas- Fermi approxima- 
tion 1241 



e HFB (n, x) sa k\J n 2 + 277J0 ps K V 2 7^o (15) 

for small quasimomenta. The fact that it is linear clearly 
shows us that, in this approximation, the low end of the 
energy spectrum corresponds to phononlike quasiparti- 
cles. Now, if we introduce ( |l4| ) into the noncondensate 
density (|l0|), we can see that the argument of the loga- 
rithm vanishes at all temperatures, making the density 
diverge at every point in space. The first line of (HG) 
shows that the divergence in the integral is caused at 
its lower limit; this restates the conclusion arrived at in 
Ref . (29) : low-energy phonons destabilize the condensate 
in the two-dimensional thermodynamic limit when the 
noncondensate quasiparticles obey the Bogoliubov spec- 
trum. 

In the Hartree-Fock approximation, on the other hand, 
the energy spectrum tends in this limit to e n 2 + 7710 
and predicts single-particle excitations whose minimum 
energy is 7^0; this can be interpreted equivalently by as- 
signing a minimum value k 2 = jng for the excitation 
quasimomentum. This cutoff is consistent with those 
proposed in the past ^ |37], [58[ and removes the 
infrared singularity in the HFB equations. 

This momentum cutoff is robust enough that it en- 
ables one to carry out Hartree-Fock calculations even in 
the thermodynamic limit: in fact, as was first found in 
Ref. |2^] , the introduction of this limit actually simplifies 
the calculations, and it is possible to find self-consistent 
solutions by simultaneously treating the noncondensate 
in the Hartree-Fock approximation and the condensate 
in the Thomas- Fermi limit 1 35 , 36 1. This model can- 
not provide realistic density profiles at every point in 
space, since the condensate density ( pd] ) has a discontin- 
uous derivative at its edge, but predicts quite reasonable 
results outside of this region, as can be seen in Fig. [I]. 



C. The Free Energy 

We have seen that in 2D the mean-field-theory BEC 
equations admit solutions both with and without a con- 
densate. The unphysical solution should be that with the 
highest free energy, since equilibrium at finite tempera- 
tures occurs when the grand potential attains a mini- 
mum; in fact, the Bose-Einstein distribution factor in 
Eq. ([?]) comes from minimizing this quantity [ pT| E5[ , 
which in our dimensionless units adopts the form 



77J0 = (/1 — x 2 — 2-yh')<d(ft, — x 2 — 2jh'), 



( 14 ) n = {U-tiN-TS)/^hcu = J d 2 x(T 



p,n-m). (16) 



where &(x) is the Heaviside step function, introduced 
to ensure that the density profile is everywhere real and 
positive. Also, since the thermodynamic limit requires 
that to — > 0, the WKB approximation becomes rigorous 
and can be used with confidence. Thus we can insert 
expression ( p"4| ) in the Bogoliubov energy spectrum (|^) 



At this point we have to keep in mind that the grand- 
canonical free energy is a function of jl, not of N; in 
order to make a meaningful comparison of these energies 
at the same TV, then, we have to compare the Helmholtz 
free energies, given by A = fl + fiN. The expressions 
given below have all been derived in the grand-canonical 
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ensemble and thus contain the chemical potential; we 
will eliminate this dependence on u by adding (IN to the 
expressions we obtain. 

In the Hartree-Fock approximation, the grand- 
canonical energy density of the system is given by plj 



T- / in = $(A+^n )$ + ^ 



?/t - 1 



-in' 2 ; 



(17) 



the first term corresponds to the condensate energy while 
the second one is the sum, weighted by the Bose-Einstein 
distribution, of the energies of the excited states; this 
last expression includes an extra term jh' 2 that has to 
be subtracted explicitly, as is usually the case in Hartree- 
Fock calculations [Q. Wc can simplify (|l7|) further by 
invoking Eq. (|J): 



T 



jJLTl 



-7n 2 



7^2 



1 



(27T) 2 



d 2 K € 
3 e 7« - 1 ' 



(18) 



Equations ( pl| ) and (]23) , of which more-general versions 
are derived in Ref. |47j by a different method fTi]] , can 
be easily seen to become identical at temperatures high 
enough that the condensate density can be neglected. 
On the other hand, their low-temperature limits differ, 
since Eq. ( p3| ) tends to a higher value than that attained 
by (pi]). This can be traced back to the fact that the 
grand-canonical ensemble has to be changed in order for 
it to correctly describe theparticle-number fluctuations 
at low temperatures [^7], Now, a 2D Bose system 

has a condensate at least at T = 0, so the two free en- 
ergies should coincide there; however, we must keep in 
mind that the semiclassical approximation used to de- 
rive Eq. ( |23| ) requires that k B T ^> Klo; thus the method 
we are using cannot describe the appearance of the zero- 
temperature condensate in the uncondensed solution. At 
the end of the next section we will study further conse- 
quences of this distinction. 



The entropy of the system can be found from the com- 
binatorial expression Q, Q 

S = -k B J2(fi log/, - {h + 1) log(/i + 1)), (19) 



which in the WKB and Hartree-Fock approximations 
yields the entropy density [El] 



1 



d 2 ne 



t 



(2tt) 2 J e £ 7* - 1 (2tt) 2 



d\ log(l-e e7 *). (20) 



The first term in (|20j) cancels with the last one in (18) 
and the other term can be integrated in closed form 
yielding 



A c = jxN - / d 2 x (j(n 2 - -h 2 ) 



-in 



92 (e 



-(x 2 +2jn-p,)/t 



)) (21) 



for the free energy of the condensed solution. (Here 
we have introduced the Bose-Einstein integral g a {x) = 
^ZfeLx x k /k". Note that gi{x) — — log(l — x), a func- 
tion that repeatedly appears above.) As T — ► 0, this 
expression reduces to the correct value in the homoge- 
neous case Jl8[ . In the thermodynamic limit, we can find 
an expression for the free energy per particle at T = 
using the Thomas-Fermi density profile: 



^TF _ 2 2^ rr— 

N 3 v 7r 7 ' 



(22) 



The free energy of the uncondensed solution is found by 
simply crossing out the condensate density in (J2l|) : 



A u =fiN 



d 2 x {^h 2 + ^-g 2 (e-^+ 2 ^ 



-A)/* 



(23) 



III. NUMERICAL METHODS AND RESULTS 

There is a time-honored prescription j2l], Q for find- 
ing the self-consistent solution of the Hartree-Fock equa- 
tions in the presence of a condensate: Initially, we assume 
that only the condensate is present and solve the Gross- 
Pitaevskii equation (||) for h' = 0. The wavefunction 
and eigenvalue that result are fed into the Hartree-Fock 
expression for the density (^), which, when integrated 
over all space, yields also a value for the noncondensate 
fraction N'; one can then readjust the condensate frac- 
tion and solve the Gross-Pitaevskii equation that results. 
The process is then iterated until the chemical potential 
and the particle fractions stop changing. 

By far the most difficult part of this process is the so- 
lution of the nonlinear eigenvalue problem (||). Different 
methods exist in the literature; we have obtained identi- 
cal results by solving it as an initial value problem |igf| 
and, much more efficiently, by employing the method of 
spline minimization |22] , p0| - This method uses the fact 
that Eq. (j^) is the Eulcr-Lagrange equation that mini- 
mizes the functional 



m = 



d 2 x [(V$) 2 + $(x 2 



-2<7n(x))<l + ^<l 4 ]. (24) 



After setting a small, nonuniform grid of fixed abscissas 
that represent the coordinate x, we take the correspond- 
ing ordinates, which represent as the parameters to 
be varied until (|24| ) attains its smallest possible value. 
Since we only have information about the value of the 
function at a discrete set of points, in order to calculate 
the necessary derivatives and to integrate we perform a 
cubic-spline interpolation; the integral is found using ten- 
point Gauss-Legendre quadrature. The minimization is 
carried out using the Nelder-Mead method. 

We checked our code by comparing its predictions in 
three dimensions to previously published results. For 
the case studied in Ref. Si . the condensate fractions we 
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FIG. 2: Total density profiles of a two-dimensional gas with 
the same parameters as in Fig. [j]. In this case we exhibit 
both the uncondensed (full line) and the condensed (dashed 
line) solutions; the latter we have broken once again into its 
condensate and noncondensate parts (dotted lines). 



found differed from those in the paper by less than one 
part in 10 . The code also reproduced previously known 
results for the ideal gas, including finite-size effects §. 

Figure [j], already discussed above, shows one of the 
solutions that we have found using the Hartree-Fock ap- 
proximation. The gas has N — 10 4 atoms; the coupling 
constant 7 has been chosen so that the system has ap- 
proximately the same radius as the three-dimensional gas 
studied in Ref. |22| , where parameters resembling those of 
the original JILA trap Q are used. The system is shown 
at T = 0.7T C , where T c is the condensation tempera- 
ture for the ideal gas. We have also found solutions for 
the Goldman-Silvera-Leggett model, which corresponds 
to the large- N limit of the GP equation. This was done 
by treating the problem as a simultaneous system of non- 
linear equations on a uniform grid and solving it with a 
least-squares method. 

When N = 10 4 , as in Fig. [J, it is not possible to find 
a self-consistent solution for the condensed equations be- 
yond T sw 0.8T C , since e HF becomes negative; this had 
already been noted in Ref. |^2| for three dimensions, 
where it was interpreted as a finite-size effect, and oc- 
curs at even lower temperatures for 2D. The limitation 
becomes more severe as N increases: for N = 10 6 we 
cannot find solutions beyond T « 0.5T C . The authors 
of Ref. HI] report predictions at temperatures very close 
to the transition by using a finite-size correction to the 
chemical potential |5l| . Our inability to work above cer- 
tain temperatures might be a consequence of using the 
semiclassical approximation, though we point out that 
the Popov approximation is expected to break down close 
to the transition temperature (TJ. [2]J . 

We also found self-consistent three-dimensional solu- 



tions using the more general Bogoliubov spectrum by 
applying the same method and using the Hartree-Fock 
solutions as a starting point for the iteration. We find 
that these solutions exhibit an enhanced depletion of the 
condensate, in agreement with those found by other au- 
thors [jsif . It was impossible, however, to find this kind 
of solution in two dimensions, even for systems with ./V as 
low as 100: after a few iterations, the chemical potential 
became too large, the Bogoliubov energies became imag- 
inary, and the noncondensate density diverged. When 
the possibility of phononlikc excitations is allowed, then, 
the condensate is destabilized, just as had been found in 
the Thomas-Fermi limit |2jj] . 

The uncondensed case, on the other hand, gives us 
solutions in these conditions, and to it we now turn. 
Equation ( |TT| ) along with the particle-conservation con- 
dition (||) for all temperatures is most easily solved by 
rewriting (11) as 



Ze~ x2/t = 2e~( 7r -^ x > sinli7rz/(x) 



(25) 



where have introduced the fugacity Z = e^l 1 and f(x) = 
2n(x) ft. Given Z, t, and 7 it is possible to find v at every 
point using a standard root-finding algorithm. One then 
wants to find the value of Z such that the total number 
of particles is N. It is better, however, to write Eq. (|2E 
at the origin, 



2e 



-(* 



^ )v ° shuWo, 



(26) 



eliminate Z between ( |25| ) and (|26|), and solve for v$ = 
j/(0), the density at the center of the trap, using the same 
root finder. 

In Fig. H we show the (Hartree-Fock) condensed and 
uncondensed solutions that we have obtained. They are 
similar in shape and exhibit identical behavior for large x. 
The uncondensed solution has a lower value at the origin 
and predicts a wider radial density profile. 

We also calculated the free energy corresponding to 
each solution; the results are shown in Fig. E| which shows 
A/N as a function of temperature for both cases when 
N = 10 3 and N — 10 4 . The free energies appear to coin- 
cide at high temperatures; this was to be expected, since 
Eqs. (|2l| ) and ( p3| ) become identical at temperatures high 
enough for the condensate density no to be neglected. As 
for the low-temperature limit, we have already noted that 
the free energy of the condensed solution tends to the 
value predicted by the zero-temperature Thomas-Fermi 
limit, while that of the uncondensed solution tends to a 
higher value. This value can actually be calculated: it is 
easy to show [pT| that the low-T limit of Eq. (pa) is 



27n(x) = /I — x 2 , 



(27) 



the Thomas-Fermi limit but with 7 replaced by 27; a 
larger interaction strength, as we can see from (|22|), im- 
plies a higher free energy. Interestingly, when we com- 
pare a given Hartree-Fock solution to an uncondensed 
solution with half the interaction strength, we find that 
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the density profiles become very similar (and are indis- 
tinguishable at T — 0), while the free energies coincide 
almost exactly for a wide range of temperatures; attrac- 
tive as this possibility might be, however, it has a serious 
flaw: the free energies start to differ as the temperature 
increases, where they should coincide by definition. This 
tells us that the factor of 2, which turns out to be the 
same one discussed after Eqs. (|]) and (p3|), has to be 
retained; the WKB approximate expression for the un- 
condcnsed state is valid only for k B T ^> Tiuj, so dropping 
the factor of 2 in order to match the T = conden- 
sate is invalid. (The authors of Ref. J3lj , in fact, omit 
this factor from their paper, although they address this 
question in a subsequent publication |32|). It is in fact 
this factor that guarantees that the uncondensed solu- 
tion has a higher free energy than the condensed one at 
all temperatures; this, despite our inability to find con- 
densed solutions using the Bogoliubov energy spectrum, 
leads us to conclude that, at least at this level of approx- 
imation, the uncondensed solution is unphysical and the 
two-dimensional finite trapped system will exhibit some 
sort of condensation at finite temperature. 

IV. CONCLUSION 

We have found solutions to the two-dimensional HFB 
equations in the Hartree-Fock approximation, both for 
finite numbers of atoms and in the thermodynamic limit; 
still, when we try to go beyond this scheme and con- 
sider the whole Bogoliubov excitation spectrum, the low 
end of which is described by phonons, we are unable to 
find self-consistent solutions for finite — even low — values 
of TV. We have seen that it is possible to describe the sys- 
tem as having no condensate at all, but these solutions, at 
least for the parameter combinations that we have stud- 
ied, have a higher free energy than their partly-condensed 
counterparts; this leads us to conclude that the 2D sys- 
tem will have some kind of condensate at low enough 
temperatures, insofar as the Hartree-Fock approach can 
successfully describe the quasicondensate that we expect 
to find in a finite system. 

One could argue that the whole mean-field approach 
we have adopted is wrong in two dimensions, since a con- 
densation into a single state is being assumed from the 
start. However, an alternative treatment |5^] that does 
not make this assumption ends up with equations iden- 
tical to (|J), so we are left none the wiser. We have also 
bypassed the fact that the interaction strength g is not re- 
ally constant in 2D, but rather depends logarithmically 
on the relative momentum |n], p3[ ; this, however, has 
been found to be of little consequence |32| . Another pos- 
sibility is that condensation occurs into a band of states, 
forming a "smeared" or generalized condensate J54|, B5| ; 
this alternative situation is not ruled out by the standard 
proof of Hohenberg's theorem |36| . 

Monte Carlo simulations do seem to predict the pres- 
ence of a condensate in two dimensions [1331 |33. Fur- 
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T/T 
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FIG. 3: Free energies per particle for the condensed (dashed 
line) and uncondensed (full line) solutions. The curves with 
higher value at T — correspond to N = 10 4 atoms, while 
those with the lower value at T = correspond to N = 10 4 . 
The coupling constant is 7 = 0.1 in both cases. The open 
circles on the vertical axis are the Thomas-Fermi predictions 
given by Eq. (jzj) for T = 0. The free energies seem to be con- 
verging in the high-temperature limit; the low-temperature 
limit, on the other hand, shows the discrepancy that results 
from restricting condensate fluctuations (see discussion af- 
ter Eq. ([23])). Note that T c , the transition temperature for 
the ideal 2D trapped gas, depends on TV; thus the x axis for 
each TV represents a different actual temperature. 

thermore, two-dimensional condensates appear to have 
been produced in the laboratory [Q Presumably 
the MC simulations show both the effects of a quasi- 
condensation of a finite system and the BKT transition 
to the superfluid state. Our attempt here has been to 
test the possibility of representing these Bose effects by a 
relatively simple set of HFB or Hartree-Fock equations. 
While no solutions can be found for the HFB set, even 
for a finite system, the Hartree-Fock equations do pro- 
vide a description of a condensed state. We feel there is 
reason to believe that this description should be a fair 
representation of the actual situation. 
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